Hanbury Brown-Twiss interferometry with electrons: Coulomb vs. quantum 

statistics* 
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A longstanding goal of Akira Tonomura was to observe Hanbury Brown-Twiss anti-correlations 
between electrons in a field-emission free electron beam. The experimental results were reported 
in his 2011 paper with Tetsuji Kodama and Nobuyuki Osakabe [l|. An open issue in such a mea- 
surement is whether the observed anti-correlations arise from quantum statistics, or are simply 
produced by Coulomb repulsion between electrons. In this paper we describe a simple classical 
model of Coulomb effects to estimate their effects in electron beam interferometry experiments, and 
conclude that the experiment did indeed observe quantum correlations in the electron arrrival times. 

PACS numbers: 



I. INTRODUCTION 

Since the pioneering detection by Hanbury Brown and 
Twiss (HBT) of "bunching" of photons in a light beam 
0, HBT experiments with massive bosons, e.g., atomic 
beams 0, [3 and ir and K mesons in high energy nuclear 
collisions [5, H[ have shown similar two-particle correla- 
tions. Seeing anti-correlation - or anti-bunching - effects 
in experiments with identical fermions where the two- 
particle intensity (/) correlation function 



C(r 2 - ri,t 2 - *i) = 



(1) 



should fall, at small separations (either in position or mo- 
mentum space) , to zero for particles of the same spin (or 
to 1/2 for unpolarized particles) has proven more elu- 
sive. Experiments with neutrons are complicated by a 
low energy nuclear resonance, while experiments with 
protons are complicated in addition by Coulomb repul- 
sion @, H- O n the other hand, anti-bunching of neu- 
tral cold fermionic 40 K atoms |{|, and the correspond- 
ing bunching of neutral cold bosonic 87 Rb atoms [T3| 
emitted from optical lattices has been successfully ob- 
served. In adddition, anti-bunching with neutral atomic 
3 He beams, as well as bunching with neutral atomic He 
atomic beams, was clearly demonstrated in the experi- 
ment of Jeltes et al. Q . 

Detecting anti-bunching in a beam of electrons has 
been a major experimental challenge over the years, ow- 
ing to the low degeneracy as well as the short coherence 
time of the beams. Starting in the 1990's Akira Tono- 
mura and his group focussed on seeing this striking effect 
of quantum statistics with electrons in a field-emission 
electron beam. Following his group's theoretical feasi- 
bility analysis [ll], E3 s electron HBT experiments have 
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been realized in free space 0, [l3j ; such experiments with 
electrons show a reduction in the correlation function for 
small space-time separation, generally attributed to anti- 
bunching. On the other hand, repulsive Coulomb inter- 
actions between electrons also reduce the probability of 
two electrons being close in space. Whether the observed 
anti-bunching effect is due to electron quantum statistics 
or rather Coulomb repulsion is the issue we deal with in 
this paper. We conclude that the recently reported ex- 
periment of Kodama, Osakabe, and Tonomura [l[ very 
cleanly sees HBT in the arrival time correlations of elec- 
trons pairs; In the experiment of Ref. 14, at a significantly 
lower beam energy, Coulomb effects account for several 
percent of the HBT signal. 

HBT interferometry has in fact been seen for electrons 
in semiconductor devices, [rH[l5j where, owing to screen- 
ing, Coulomb effects arc less important. For example, in 
the HBT experiment of Ref. 16, the screening length 
(~ 5nm) is typically much smaller than the Fermi wave- 
length (~ 40nm) [lq . Interestingly, two dimensional 
mesoscopic semiconductors open the possibility of see- 
ing HBT correlations for fractional statistics [T3, [H| as 
well as Aharonov-Bohm physics [l9[ . 

The importance of correcting for Coulomb interactions 
has long been reco gniz ed in interpreting high energy nu- 
clear experiments [20| . For example, the raw correla- 
tion function of distinguishable pions of opposite charge 
(tt + it~ ), produced in the E877 ultra-relativistic heavy ion 
collision experiment 21| , shows a very similar buildup at 
small momentum differences to those of identical charged 
pions (7r + 7r + and 7r~7r~). The Coulomb interaction be- 
tween opposite charges tends to increase the probability 
of a pair of bosons being close in momentum, while re- 
ducing that for like charges. Only after the effects of 
Coulomb interactions are extracted, does one see the ex- 
pected effects of quantum statistics (see, e.g., Refs. 7 and 
22). 

Our aim in this note is to present a simple schematic 
discussion of Coulomb effects in interferometry experi- 
ments with electrons, based on the classical behavior of 
electrons taken pairwise. We do not attempt to explain 
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FIG. 1: Schematic of emission of two electrons. The first 
(closed circle) and second (open circle) travel downward to 
the detector plate at velocity ~ v. Coulomb acts along the £ 
direction. 

the detailed results of Tonomura's group on HBT with 
electrons, but rather aim to estimate the role of Coulomb 
interactions in their search for quantum correlations. 

The Coulomb problem for a pair of electrons is charac- 
terized by four length scales, 1) the electron Bohr radius 
ao = h 2 /me 2 , with m the electron mass; 2) the size r$ of 
the emitting region transverse to the beam; 3) the typical 
separation zq of the particles along the beam direction; 
and importantly, 4) the classical turning point of the pair, 
r^, defined by e 2 /rt p = q 2 /2m re d, where q is magnitude 
of the final relative momentum of the two particles and 
m red — m/2. 

The traditional method of correcting for Coulomb in- 
teractions is to employ the Gamow correction, which 
assumes that the characteristic separation of the pair 
of particles is much smaller than their classical turning 
point, namely, that the particles are produced well within 
the classically forbidden regime bounded by r tp . The ac- 
tual rate observed in an experiment is taken to be that 
in the absence of Coulomb interactions times the Gamow 
correction, |?/> c (0)| 2 , which is the absolute value squared 
of the relative Coulomb wave function at the origin, 

M»=(^)" 2 . (2) 

where the dimcnsionless parameter r\ equals zz'e 2 /hv re i, 
with Vrei the relative velocity of the two particles of 
charges ze and z'e. On the other hand, in a field emission 
source, the relative energy of a pair is ~ AE, where AE 
is the initial electron energy spread in the beam. Thus 
r tp ~ e 2 /AE ~ (1.46/ AE eV ) nm, where AE eV is the 
energy spread measured in electron volts. Since r tp is 



FIG. 2: Schematic picture of the effect of Coulomb repulsion 
on the trajectories of two electrons emitted successively from 
the same point in the tip. 



generally tens of nm, r tp <C Tq, indicating that a pair 
of electrons in such an experiment is typically emitted 
outside the pair's classical turning point. Coulomb ef- 
fects are dominated by classical physics rather than by a 
quantum Gamow correction. 



II. CLASSICAL MODEL 

To bring out the effect of Coulomb interactions, we 
model the experiment as independent emissions of elec- 
trons from a tip, followed by acceleration to final velocity 
v in the beam direction (z) and energy Ef = mv 2 /2. 
We first neglect quantum statistics, and focus on the 
Coulomb effects in a single pair of particles, since the ma- 
jor contribution to the correlation function is from par- 
ticles nearby in space and time, a configuration in which 
we can, to a first approximation, neglect many-body ef- 
fects. We assume that the emission points of the pair are 
separated by xi in space and tj in time. Thus the initial 

1/2 

spatial separation of the pair is = (x 2 + (vti) 2 ) . 

The Coulomb repulsion between the electrons increases 
their relative separation in space and time, as is illus- 
trated schematically in Fig. [5] for two electrons emitted 
at the same point in the tip at times tA and then ts- 
The electron separation at later time is readily calculated 
from the conservation of energy of the relative motion of 
the two electrons, E re i — ms 2 /4 + e 2 /s. A lower bound 
on the size of the Coulomb hole can derived by neglecting 
the initial relative kinetic energy of the pair; then after 
integration of s, one finds that the final separation of the 
pair, Sf, is determined implicitly by 

A* = s- /2 {v / ^ r T) + hi [y^ + v^l]}, 

(3) 

where a = s//sj, and At ~ L/v is the time elapsed be- 
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twccn emission and detection of the pair. Approximately, 

{Sc (^c/^i) i &i ^ $c ( a \ 

(4) 
Si-) Si Sc; 

1 /3 

where s c = vr c = (2e 2 L 2 /Ef} defines a critical 
Coulomb distance and time r c ; numerically, s c ~ 6.5 x 

10- 4 i^ 3 /^}5ey cm - 

As seen in the left panel of Fig. [3j the relation be- 
tween the final separation Sf when the electrons reach 
the detector plate and Sj is not monotonic, but rather is 
decreasing at small Sj and increasing at large s, . For very 
small initial separation, the Coulomb interaction signifi- 
cantly accelerates the two electrons away from each other, 
making the final separation large, while for very large ini- 
tial separation, the Coulomb interaction is negligible, and 
the final separation is essentially equal to that initially. 
No matter how small is, the final spatial separation 
is finite, i.e., there is a Coulomb hole in the distribution 
of final separations. Coulomb forces increase the spatial 
separation between a pair of electrons, s(t), with time, 
so that Sf > Si. 

Since the angle 9 that the relative position vector of 
the electrons makes with respect to the beam axis (see 
Fig. [T]) is conserved in the motion, the final separations 
at detection are related by Xf/xt = i//tj = s//sj; thus 
the final separation in time of the two electrons is given 
by tf — (sf/si)ti. At small ti/r c , the final tf has the 
structure shown in the right of Fig. |3l dependent on Xi. 
However, for ti/r c ^ 0.1, the curves on the right of Fig. [3] 
converge to that in the left panel of Fig. [31 which has 
a minimum at Sf f=s s c . To a first approximation, the 
minimum i/ ~ s c /v = t c is independent of the initial 
spatial separation. 

Experimentally one measures the correlation function 
C(tf) in terms of the subsequent time intervals tf be- 
tween detection of particles, averaged over the distribu- 
tion of initial emission intervals ti . For independent emis- 
sions, the distribution of times between adjacent emis- 
sions from the tip is Poissonian 

P (U) = ie-Wfc ( 5 ) 

where ti is the average time separation between two emis- 
sions. In the absence of Coulomb corrections and quan- 
tum statistics, Co(ii) = 1. The final C{tf) and P{tf) are 
given in terms of the map ([3]) between s/ and s^ since 
the map is not simply one-to-one, one needs to sum over 
the two branches. The resulting P{tf) and C(tf), cal- 
culated with the approximate solution (T?|) are shown as 
thin lines in Fig. @] 

In experiments, the finite time resolution of the de- 
tectors would smooth out the sharp Coulomb holes in 
Fig. |4] Measurement of an observable f(t) at time t av- 
erages the actual f(t') over t', weighted by the detector 
time resolution function R(t — t'): 

/+oo 
R(t-t')f(t')dt', (6) 
-OO 
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FIG. 3: a) Final vs. initial spatial separations, with Coulomb 
interactions included classically, and b) final vs. initial time 
separations for different initial transverse separations Xi = 0, 
5, 10, and 30 nm (top to bottom). For large ti/r c , the curves 
all converge to that in the left panel. 
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with f(t) = f(-t) for t < 0. In Fig. GO effects of time 
resolution are indicated by the thick curves, where we 
take a Gaussian time resolution function 

R { t-t') = ^-e-^') 2 ^ , (7) 

with a characteristic time scale t r chosen here for illustra- 
tion to equal t c . Typically, t r 3> t c . The dip at low t f in 
right panel illustrates clearly how Coulomb correlations 
can mimic quantum correlations. 



III. HBT ANTI-CORRELATIONS 

To assess the importance of the Coulomb hole it is nec- 
essary to compare it with the regime of suppression of 
the correlation function from quantum statistics. The 
Pauli principle suppresses the correlation function be- 
tween same spin electrons emitted from nearby points 
on a time scale 

tHBT = H/Tf, (8) 
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FIG. 4: Effects of Coulomb repulsion as seen in a) the final 
time distribution P(tf) vs. tf, with U = 0.2 ns; the vertical 
axis is in units of 1/U, and b) the normalized second-order 
final correlation function C(tf) vs. tf. The thin lines are 
calculated directly from Eq. Q; the thick lines include finite 
time resoluton via Eqs. ([6]) and {?]), with t r = r c . 



where Tf is the effective longitudinal temperature of the 
electron gas after acceleration. Since the acceleration of 
the beam is essentially adiabatic, the entropy per electron 
is conserved, and the temperature of the gas falls with 
expansion. Owing to acceleration the density n of the gas 
drops, since in a steady state the current nv remains con- 
stant from emission through acceleration. To estimate 
the final gas temperature, we note that the entropy per 
particle of a gas with an anisotropic temperature depends 
on T x / 2 Tx/n, where T is the longitudinal and T± the 
transverse temperature. Thus in anisotropic expansion, 
T 1 / 2 T±/n remains constant, and for T± constant, Tv 2 is 
invariant. Initially v = \jTijm and after acceleration to 
velocity much greater than that of the thermal motion, 
v = W2Ef/m. Thus the final longitudinal temperature 
of the beam is given by Tf y/2Ef/m ~ T y/Ti/m, so that 
with Eq. © and T 4 ~ 2AE, we have T f = 2(AE) 2 /E f 
(essentially the result derived heuristically in Ref. 1) and 
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(9) 



The ratio of the Coulomb to HBT suppressions is thus 

,11/6 



tHBT 
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7 2 r 2 / 3 ' 

^eV cm 



(10) 



This simple calculation indicates the importance, in an 
HBT measurement of correlations in arrival times, of ac- 
celerating the electrons to a large final beam energy in 
order to overcome Coulomb effects. 

In addition to time anti-correlations in the beam, HBT 
correlations should appear in the electron spatial separa- 
tions, analogous to the spatial correlations seen in the 
original Hanbury Brown-Twiss measurement of the an- 
gular diameter of the star Sirius @] . When two electrons 
are emitted at the same time, but spatially separated, 
the correlation function is suppressed in space on a scale 
of order the particle wavelength divided by the angular 
size of the source, or 
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f,KeV r Wnm 



where rio nm is the transverse size of the emission region 
in units of 10 nm. Comparing with the minimal Coulomb 
hole, s c , we find 



srbt 



L 



1/3 



E 



1/6 

f,KeV 



(12) 



r 10r, 



To reach this ratio, the initial spatial separation must 
be of order s c ; however, a more realistic bound on s 4 is 
the transverse size of the emission tip, which leads to a 
considerably larger Coulomb hole, as one can infer from 
Fig. 1. 

In the experiment of Kodama, Osakabe, and Tono- 
mura Q E f ~ 50-100 KeV, and AE ~ 0.17 cV; with 
L ~ 100 cm one estimates that tHBT /t c ~ (2 — 6) x 10 3 , 
sufficiently high that one does not need to worry about 
Coulomb effects in measuring pure time-of-arrival corre- 
lations. On the other hand, in the opposite regime, when 
measuring spatial HBT correlations, one has shbt/s c ~ 
0.5, indicating that Coulomb effects must be taken into 
account in analyzing the experiment. 

In contrast, in the lower energy experiment of Kiesel 
et al. [HI, where E f - 0.9 KeV, AE ~ 0.13 cV, and L ~ 
1 cm, one has tHBT /tc ~ 44, and thus Coulomb effects 
while small are not entirely negligible. For spatial cor- 
relations, however, shbt/sc ~ 0.25. and thus Coulomb 
effects are dominant. 

In conclusion, a full analysis of the HBT experiment of 
Kodama, Osakabe, and Tonomura requires correcting for 
Coulomb effects. The simple model presented here forms 
a useful and simple basis for including Coulomb interac- 
tions among the electrons. A detailed analysis requires 
taking into account the distribution of initial spatial sep- 
arations and velocities of the electrons in addition to the 
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time separations, the effects of realistic finite time reso- 
lution, as well as the effects of the quadrupole magnets 
which give one the freedom to adjust the angular size 
of the beam. While Coulomb effects are present inde- 
pendent of the relative spin of the electron pair, Pauli 
quantum correlations occur only between same spin elec- 
trons; thus to distinguish optimally Coulomb repulsions 
from quantum correlations one would ideally like to re- 
peat the experiments with spin polarized electron beams. 
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